Low energy spin wave excitations in bilayered magnetic manganite La 2 _ 22 ;Sri +2x Mn207 

(0.30 < x < 0.50) 

H. MartinhcQ and C. Rettori 
Institute) de Fisica "Gleb Wataghm", UNICAMP, 
13083-970, Campinas, SP, Brazil. 



D. L. Huber 

University of Wisconsin, Dept. Phys., 
Madison, Wisconsin, 53706, USA. 



J. F. Mitchell 

Materials Science Division, Argonne National Lab., Argonne, Illinois, 60439, USA 



S. B. Oseroff 

San Diego State University, 
San Diego, California, 92182, USA. 



We have studied the low-temperature behavior of the magnetization and the specific heat of 
the bilayered perovskite system La2-2xSr\+2 X Mn207 , for 0.30 < x < 0.50. Our analysis reveals 
that below 30 K the temperature dependence of the magnetization of the ferromagnetic samples, 
x = 0.30, 0.32, and 0.36, in a field of IT can be interpreted in terms of the thermal excitations of a 
two-dimensional gas of ferromagnetic magnons. The specific heat in zero field for these samples as 
for the x — 0.50 antiferromagnetic one, is linear with temperature between the range of 1.8K < T 
< 10K. This behavior can be also explained by the magnon gas model. By comparing specific heat 
measurements in zero field with those taken in a field of 9T we are able to extract the lattice and 
electronic contributions and determine the in-plane exchange interactions. That are found to be 
in reasonable agreement with the values inferred from the analysis of the magnetization data and 
also with the values reported by inelastic neutron scattering studies. In addition, we found that the 
electronic density of states obtained for the x = 0.50 sample is in agreement with previous band 
structure calculations. 

PACS numbers: 75.47.De;75.30.Ds;75.30.Et;75.40.Cx 



I. INTRODUCTION 



Manganese oxides with perovskite structure have been 
subject of intense study in recent years. This is in large 
part motivated by two main factors: their potential tech- 
nological application, related to the Colossal Magnetore- 
sistance (CMR) effect displayed by many of these oxides, 
and the large variety of magnetic/electronic phenomena 
such as double-exchange, superexchange, Jahn- Teller ef- 
fect, charge/orbital ordering related to the CMR. 0,0 

A common way to represent the manganese per- 
ovskite oxides is through the Ruddlesdcn-Popper se- 
ries (R,A)„ + iMn„03„+i (R = rare-earth; A = alkaline 
metal). The basic structure of the series consists of al- 
ternate stacking of rock-salt type block layers (R,A)2 02 
and n— Mn02 sheets along the c— axis, where n rep- 
resents the number of adjacent Mn02 sheets, a quan- 
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tity closely related to the effective dimensionality of the 
system. So, the n — oo series represents the proto- 
type three-dimensional (3D) CMR material (R,A)Mn03, 
the subject of many investigations lately; n — 1 corre- 
sponds to compounds such as La2Cu04, the metal ox- 
ide that originates the high-Tc superconducting cuprates 
family; n = 2 represents the bilayered magnetic man- 
ganite system (R,A)aMn207. This last system shows 
two-dimensional (2D) magnetic and electronic proper- 
ties associated with a very large magnetoresistance (MR) 
p(0) /p{H) ~ 20, 000% at T c ~ 130K and H = 7T for 
Lai^Sr! gMn 2 07.:3i| For comparison, in the 3D Sr-doped 
La .85Sr . 15 MnO3 the MR reaches ~ 110% at H = 15T 
around T c .[| 

The structure of the n = 2 bilayer system is shown in 
figure^ It consists of pairs of corner shared MnO@ octa- 
hedron stacked along the c— axis and isolated by La(Sr) 
atoms. The dashed lines enclose the unit cell of the 
bilayer system. Also, indicated in the figure are the 
intra-Mn02 plane, Jm, inter-Mn02 planes, Jj_, and inter- 
bilayers, J', exchange constants. 
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FIG. 1: Structure of the hole doped bilayered system 
La2-2xSri+2a;Mn207. Dashed lines represents the unit cell 
of the bilayer where J\\ , J± and J' are the exchange constants 
described in the text. 



La2-2a;Sri + 2 2: Mn207 present a wide variety of elec- 
tron transport behaviour and magnetic structures. For 
the range 0.30 < x < 0.50 it is found at low tempera- 
tures, T < 10K: (i) metallic character and ferromagnetic 
(FM) order along the c— axis for x — 0.30; 5] (ii) metallic 
character and FM order in the MnC>2 sheets (ab— plane) 
for x = 0.32; [j| (hi) metallic character and FM order in 
the ab— plane, with a weak antiferromagnetic (WAFM) 
component along the c-axis for 0.32 < x < 0.48; |5j (iv) 
metallic A-type AFM with a small fraction (< 18%) 
of a secondary CE-type AFM insulating phase for 
0.48 < x < 0.50.0 

In all cases, the 2D character of the magnetic interac- 
tions is evident in these systems. For example, the mea- 
sured critical exponent of the sublattice magnetization 
near Tc in the x = 0.40 sample, was found to be very 
close to the expected value for the 2D Ising model. [8j The 
two-dimensionality of these systems is also revealed by 
the strong anisotropy in the resistivity, with the ab— plane 
resistivity being two or three orders of magnitude smaller 
than that along the c— axis.0 

We have studied the influence of the dimensionality on 
the magnetization and specific heat of the n = 2 hole- 
doped bilayered system La2_2i C Sri + 2: i; Mn2C , 7 (0.30 < 
x < 0.50). We present a consistent interpretation of the 
low-temperature data in terms of a gas of thermally ex- 
cited 2D magnons. 



II. EXPERIMENTAL DETAILS 

Single crystals of L^^rrSri^zM^Or for 0.30 < x < 
0.50 were melt grown in flowing 20% O2 (balance Ar) in 
a floating zone optical image furnace (NEC SC-M15HD). 
Magnetization measurements have been taken in a Quan- 
tum Design dc SQUID MPMS-5T magnetometer. The 
specific heat was measured using the specific heat insert 
of a Quantum Design PPMS-9T measurement system. 

III. EXPERIMENTAL RESULTS 

Figure [21 shows the field-dependent dc— magnetization 
curves, MvsH, taken at 10K for H || ab and H || c. 
The measurements in all samples were made by zero field 
cooling (ZFC) to 10K, then the field was increased to 
5T and decreased again to zero. For x = 0.30, 0.32 and 

0. 36 we observed the usual non-linear FM behavior. The 
x = 0.50 sample displays a linear AFM behavior with a 
slight anisotropy between the c— axis and the ab— plane 
plus a small non-linear FM-like behavior at low fields 
[H < lOkOc). 

From the MvsH curves we obtain the saturation mag- 
netization, Ms, for the FM samples. Using this values it 
is possible to estimate the effective spin of the Mn ions, 
since S e ff = Ms/2(J.b, and compare it with the theoret- 
ical expected one, Sfff. 

Considering that the theoretical valance of the 
manganese ion in the undoped compound x = 0, 
La2SrMn 2 07, is 3+, and that their electronic configu- 
ration is 3d 4 4s°, their expected average effective spin is 
S£iy(x — 0) = 2. Under hole-doping this value will be 
reduced by x/2, giving S*J}f(x) = (4- x)/2. On Table 

1, we note that the agreement between S e f / and S 1 * h y is 
at best reasonable. That indicates that the average spin 
model does not provide a very good description of the 
saturation magnetization. 



TABLE I: M s , S eff and Sfff for the FM samples. 



X 


M s (emu/g) 




Qtheo 
b e.f.f 


0.30 


68.5 


1.63 


1.85 


0.32 


77.3 


1.92 


1.84 


0.36 


75.3 


1.86 


1.82 



Figure shows the temperature dependence of magne- 
tization, MvsT curves, for the FM samples. All samples 
were field cooled from 300 to 1.8K in a IT applied field. 
The curves show a FM transition at Tc ~ 100 — 140K, as 
reported previously, |2| with a slight anisotropy between 
H || c and H \\ ab, which approaches zero for x = 0.32. 

In the left scale of the figure 0] the MvsT curves for 
the AFM sample are shown. The sample was field cooled 
from 750 to 1.8K in a IT applied field. As pointed out 
by Kubota et al.^ when analyzing their magnetic neu- 
tron diffraction data, this sample presents two kinds of 
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AFM phases, an A-type and a CE-type. The Neel tem- 
perature of the A- AFM phase coincides with the charge 
ordering temperature at Tjy co — 210K, while the CE- 
AFM orders at T% E = 145K.@ These two characteristic 
temperatures are shown in figure 01 by vertical dotted 
lines. 

The right scale of figure 0] shows the temperature de- 
pendence of the inverse of magnetic susceptibility, \ ■> 
for the x = 0.50 sample for both directions, H || c and 
H || ab. Above T ~ 100K the two curves are identical. 
The number of effective Bohr magnetons, p e ff, and the 
Curie- Weiss temperature, 9, extracted from fitting the 
data tox^ 1 = {T- 9)/C, for the range 600 < T < 750K 
are p eff = V^C = 4.59^s/Mn and 6 = 300K. For 
this compound, the expected valance of Mn is 3.5 and 
Peff — 4.4/ie, in good agreement with our experimental 
results. 

Figure |3] shows the specific heat results plotted as 
C/TvsT 2 , at zero field and H = 9T, applied parallel to 
the c— axis for all samples. Two interesting features can 
be observed. The first one is the linear dependence of the 
zero-field data that shows a high value for gamma. This 
feature is basically independent of the hole concentration, 
magnetic order and conductivity. The second feature is 
the strong magnetic field dependence of the data. In the 
next section we will show that the origin of this behavior 
results from the 2D character of the system, namely, the 
influence of the 2D magnons. We will analyze the influ- 
ence of the 2D magnons on the magnetization of the FM 
samples and on the specific heat of all compounds. 
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FIG. 2: MvsH measured at 10K for x = 0.30, 0.32, 0.36 and 
0.50 samples. Solid symbols correspond to H || ab and open 
symbols to H \\ c. 
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FIG. 3: MvsT for the FM samples x = 0.30, 0.32 and 0.36 
taken at H = IT. The solid symbols correspond to H \\ ab 
and open symbols to H || c. 



IV. THEORETICAL MODEL 

A. Magnetization 

To analyze the change on the magnetic response due 
to elementary excitations, it is useful to consider the 
fractional change of the magnetization, [M(T — > 0) — 
M(T)]/M(T -> 0) = AM/M(0). Using boson statistics, 
it can be expressed as:0 

AM _ S - (S z ) _ J_ f°° kdk 

M(0)~ S ~inSj e %&_{ ' ' 

for spin waves in a system with both an acoustic and an 
optical branch, with the latter not thermally populated. 
Here S is the spin, (S z ) is the average value of the z 
component of the spin, k is the magnon wave vector and 
huj(k) is the magnon energy. 

Considering 2D FM acoustic magnons, with wave vec- 
tor k = k a +k;, (measured in units of the reciprocal of the 
nearest-neighbor distance, 3.87A), the dispersion relation 
at low temperatures, neglecting the interbilayer exchange 
energy J', in the presence of a field, can be written as 

fiw(k) = ~J {l Sk 2 + 2[i B H = Dk 2 + A, (2) 

with D = — JuS being the spin wave stiffness constant 
and A = 2hbH the gap induced by the magnetic field. 
Evaluating using the dispersion relation we ob- 
tain 
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FIG. 4: Le/t scate: MvsT for the AFM sample x = 0.50 
measured at H = IT. The vertical dotted lines represent the 
Neel temperatures of the two AFM phases at T$ E = 145K 
for the CE type AFM structure and T^ co = 210K for the 
A type one; Right scale: x~ lvs T for the AFM sample x = 
0.50, with the solid line showing the Curie- Weiss fitting for 
the parameters p e ff ~ 4.59/^s/Mn and 8 — 300K. The solid 
symbols correspond to H \\ ab and the open symbols to H \\ 
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with A in units of K. 

The corresponding analysis for A-AFM magnons is 
more complicated and will not be discussed here. 

B. Specific Heat 

The contribution to the specific heat from the 
magnon excitations is calculated summing over harmonic 



oscillator-like contributions 



Huj_\2 
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We can 



(e*B' r -1) 

express the total contribution to the specific heat coming 
from the acoustic branch of the 2D FM magnons having 
the dispersion relation J2J as 



C n 



4irk B T 2 J 



k[hw{k)] 2 e "b t 
[e fc s T - l] 2 



dk 



(4) 



Using this expression we can derive the zero-field, A = 
0, 2D FM magnon specific heat as 
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FIG. 5: C/TvsT 2 for x = 0.30 (a), 0.32 (b), 0.36 (c) and 
0.50 (d) samples. Open symbols represent the zero field mea- 
surements. Solid symbols corresponds to measurements with 
H = 9T || c. 



where "f ma g — ~KRks/2AD. This expression was first de- 
rived by Colpa|ll| for a 2D Bravais lattice of spins. The 
field-dependent specific heat is obtained solving Q with 
A 7^ 0. Unfortunately, in this case, the result is not 
analytical. 

Neutron scattering results by Hirota et al. 12] indicate 
that the coupling between the two MnC>2 sheets, J±, is 
nearly zero for x = 0.50. Thus, in the analysis of the 
specific heat for x = 0.50, we neglect J± so that the 
magnetic contribution to the specific heat is the result 
of 2D excitations that are effectively FM for energies ~ 
ksT. We can then use the expression (jSJ) keeping in mind 
that in this case "f mag = irRks/^D, since there are two 
degenerate magnon branches for each fc-value. 

Considering the contributions to the specific heat com- 
ing from electrons, phonons, and magnons, the total zero- 
field specific heat can be described as 



(H = 0) = 7e + lmag + I3T 2 = 7e// + /3T 2 , (6) 



where we have defined an effective gamma-term J e ff = 
le+lmag- So, the zero field specific heat can be displayed 
as a linear plot in a C/TvsT 2 curve. 
For H ^ the total specific heat will be 



^(F^0)= 7e + /3T 2 - 



k[hw(k)] 2 e~ T ^ r 



**k B T* J _ 1]2 



dk 



( 5 ) with hw{k) = Dk 2 + A. 
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FIG. 6: AM/M(0) plotted as function of -Tln[l - e^T] for 
x = 0.30 (a), 0.32 (b), and 0.36 (c). The magnetic field of H = 
IT (A = I.34K) was applied along the easy-axis. The vertical 
dotted line indicates the temperature of 30K. The solid line 
represents the fitting of AM/M(0) to -bTln[l - e~T]. 



For x = 0.50, the prefactor in front of the integral in 
cq.Q need to be multiplied by a factor of 2 to account 
for the two-fold degeneracy mentioned previously. 

In the next section we will apply these results to the 
analysis of the magnetization and specific heat data. 



literature, 0, 0, Q, 0] we use the average theoreti- 
cal effective spin Sl'jy, instead of the experimental value 
S e ff. Using Sfff we obtained the stiffness constant,/) , 
and the exchange energy Jy given in Table II. No signifi- 
cant differences are found if S e fj is used instead. In that 
case the values of Jy are: -78(4), -52(2) and -74(4)K, 
for x =0.30, 0.32 and 0.36, respectively. Our values of 
Jy are in good agreement with the Jy values obtained 
from neutron diffraction measur ements[Il 13 El 13. 
Thus, from our analysis we can conclude that the main 
contribution to the temperature dependence of the mag- 
netization of the FM bilayer system comes from the 2D 
FM magnons. 



B. Specific heat 

Figure shows the zero field specific heat data fit to 
eq.@). In the temperature interval between 1.8K< T < 
10K, the data are in a good agreement with the FM 
magnon model. From this analysis we derived the Jeff 
and (3 values for all concentrations. 



TABLE III: Experimental specific heat parameters. 



X 




7e 


P 


D(H = 9T) 




(mJ/mole K 2 ) 


(mJ/mole K 2 ) 


(mJ/mole K 4 ) 


(K) 


0.30 


28(1) 


7.0(2) 


0.27(3) 


88(4) 


0.32 


25(1) 


7.0(3) 


0.31(3) 


96(2) 


0.36 


23(1) 


7.0(2) 


0.27(2) 


98(6) 


0.50 


24(1) 


4.0(8) 


0.40(4) 


196(16) 



V. ANALYSIS AND DISCUSSION 

A. Magnetization 

Figure shows the fractional change of the magneti- 
zation AM/M(0) as function of -Tln[l - e~£] for the 
FM samples. A field of H = IT (A = 1.34K) was ap- 
plied parallel to the easy-axis of the samples: || c for 
x = 0.30 and || ab for x = 0.32 and 0.36. We fit the data 
to AM/M(0) = — bTln[l — e~r] and the parameters ob- 
tained are given in Table II. 



TABLE II: Parameters 6, D and Ju , obtained as described in 
the text. 



X 


6(10- 4 i^- 1 ) 


D(K) 


J\\(K) 


0.30 


1.9(1) 


112(6) 


-60(2) 


0.32 


2.1(1) 


104(6) 


-56(2) 


0.36 


1.6(1) 


138(9) 


-76(4) 



From the expression J3J we have that b = l/8irDS. 
In order to enable us a consistent comparison between 
our model and the neutron diffraction analysis from 




x = 0.50 (d) 
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FIG. 7: Zero field specific heat data for x = 0.30 (a), 0.32 (b), 
0.36 (c) and x — 0.50 (d) samples. The solid lines represents 
the fitting of the data to eq. (6). 
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TABLE IV: Parameters obtained from the specific heat analysis as described in the text. 



X 


9{sf) 


^mag 


D(H = 0) 


Jy(# = 0) 


J|,( J ff = 9T) 




(states/ Ry f.u.) 


(mJ/mole K 2 ) 


(K) 


(K) 


(K) 


0.30 


40(1) 


21(1) 


104(2) 


-56(2) 


-47(1) 


0.32 


40(2) 


18(1) 


122(2) 


-66(4) 


-52(2) 


0.36 


40(1) 


16(1) 


136(2) 


-74(4) 


-54(2) 


0.50 


23(5) 


20(1) 


218(4) 


-124(6) 


-112(8) 



Figurc[S]shows the specific heat data, taken at H = 9T, 
fitted to eq.Q. In this fitting we used the (3 parameters 
derived from the zero field analysis. From it we obtain 
7 e and D{H = 9T) and calculated the exchange energy 
J|| (H = 9T). No appreciable difference is found in the Jy 
values when using or S e ff- For x = 0.50, it should 

be noted that the ferromagnetic approximation for the 
spin wave energies, as indicated by the linear behavior of 
the specific heat in zero field, only pertains to the modes 
that make a significant contribution in the temperature 
range of the experiment. If measurements were made at 
much lower temperatures, one would expect a temper- 
ature dependence characteristic of a weakly anisotropic 
A-AFM. 

We calculated the ^ mag from j mag = j e ff - j e . Using 
7 e and j mag we computed the electron density of states 
at Fermi level, g(ep), and D(H = 0). Finally, from 
D(H = 0) we obtain the exchange energy J\\(H = 0). 
The experimental parameters 7e//, 7e, P and D(H = 9T) 
are listed in the Table III. The values of <?(£f), 7mag, 
D(H = 0), J,, (H = 0) and Jy (H = 9T) derived from the 
data, are given in the Table IV. 

The values of j e and g(ep) of the insulating samples, 
x = 0.30, 0.32 and 0.36 are in reasonable agreement with 
those reported in the literature (see, e.g., the results of 
Okuda et al.\l]^), suggesting a half metallic behavior. In 
fact, electronic band structure calculations performed by 
Meskine et al. suggested a half metallic character for 
some bilayers. In particular, for x = 0.50 they found the 
theoretical value of g(ep) = 19 states/Ry f.u., in excellent 
agreement with our experimental value of 23 states/Ry 
f.u.. 

The values of Jy (H = 0) and Jy (H — 9T) are in rea- 
sonable agreement with the results pre viously reported 
by neutron scattering studies [12, EH, EH E3- In figure 
we compare the Jy values obtained in the present work 
and the ones reported by Hirota et aZ.E2 The J\\ (H = 0) 
values obtained from the zero field specific heat agree 
with the values inferred from the analysis of the mag- 
netization. The values of Jn (H = 9T) differ somewhat 
from J\\{H = 0). That may be a consequence of the 
presence of the double exchange mechanism since Jy is 
field-independent for a pure Heisenberg system. The 
Jy (H = 0) values are also in fair agreement with the Hi- 
rota 's ones, while the J\\(H = 9T) are in excellent agree- 
ment. The above is possible related to the differences in 
the energy scale probed by the different techniques. Our 
zero field thermodynamic measurements probe lower en- 
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FIG. 8: Fit of the specific heat data measured at H = 9T || 
c to eq.(7) for x = 0.30 (a), 0.32 (b), 0.36 (c), and 0.50 (d) 
samples. 

ergy magnons, while the field-dependent measurements 
involve magnons with a slight higher energy, as the 9T 
magnetic field opens a gap ~ lmeV in the magnon dis- 
persion relation. On the other hand, generally neutron 
experiments involve magnons with energies > 3meV. So, 
it is expected a smaller difference between the values of Jy 
derived from neutron diffraction and the obtained from 
9T specific heat than with the ones extracted from the 
zero field specific heat. Also, a significant renormaliza- 
tion of the magnon energies at small k, that is not re- 
flected in the neutron data, may be present. 

VI. CONCLUSION 

In summary, our experimental specific heat results re- 
veal two interesting features: (i) the linear, high-gamma 
behavior of the zero-field data, independent of hole con- 
centration, magnetic order or conductivity and (ii) the 
strong magnetic field dependence of the measurements. 

We propose a model considering the presence of 2D 
FM magnons to explain the data. We calculated the in- 
fluence of these magnons on the magnetization of the FM 
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FIG. 9: Comparison between the Ju values, obtained on the 
present work and the ones reported by Hirota et al.|l2j Open 
squares correspond to the values derived from the MvsT anal- 
ysis. Open (solid) circles correspond to the values derived 
from C/TvsT 2 analysis with H = 0(9T). Open triangles 
present the values from Hirota et 



samples and specific heat. In all the cases the model gives 
a good fit, showing that the change in magnetization and 
the specific heat of the system can be well described by a 
2D FM magnon gas. In the case of the antiferromagnetic 
system (x — 0.50), the linear behavior of the specific heat 
is shown to be a consequence of the weak coupling be- 
tween the bilayers sheets, that enables the occurrence of 
the 2D FM magnons. 

As noted, we find similar values for Ju from both the 
magnetization and the specific heat analysis, indicating 
an internal consistency on our analysis. Moreover, the 
values of Jii we found are close to the ones derived from 
the neutron scattering studies^, 0, 0, [^. The last 
provides further support to our analysis. 
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